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Within the non-linear Ginzburg-Landau (GL) theory, we investigate the vortex structure in a 
superconducting thin film with a ferromagnetic disk on top of it. Antivortices are stabilized in 
shells around a central core of vortices (or a giant-vortex) with size-magnetization controlled "magic 
numbers" . An equilibrium vortex phase diagram is constructed. The transition between the different 
vortex phases occurs through the creation of a vortex-antivortex pair under the magnetic disk edge. 



PACS numbers: 74.78.-w, 74.25.Op, 74.25.Dw. 

Almost half a century ago, Abrikosov used the 
Ginzburg-Landau (GL) equations to predict that the 
magnetic field penetrates a type II superconductor in 
bundles arranged in a regular lattice, i.e. the vortex lat- 
tice Q. Nowadays, it is well known that the motion of 
Abrikosov vortices gives rise to dissipation which is the 
limiting factor for the size of the critical current of a 
superconductor. The last decade has seen an increased 
interest in vortex matter in inhomogeneous superconduc- 
tors, where defects (random and ordered) are used to 
pin the Abrikosov lattice in order to increase the critical 
current which is crucial for practical applications. Irra- 
diation by heavy ions allowed to create dense random 
arrays of columnar defects which increased the critical 
current density 0] . Thanks to substantial progress in the 
preparation of magnetic microstructures in combination 
with superconductors j^j, superconductor/ferromagnet 
hybrid systems became very interesting both theoreti- 
cally and experimentally. In Ref. p] Marmorkos et al. 
investigated the problem of a magnetic cylinder with 
out-of-plane magnetization embedded in a superconduct- 
ing film. They solved the non-linear Ginzburg-Landau 
equation numerically, with appropriate boundary condi- 
tions, and found a correspondence between the value of 
the magnetization and the vorticity of the giant-vortex 
states. 

More recently, the interest shifted to the pinning be- 
havior of regular arrays of submicron ferromagnetic disks 
where additional pinning contributions arise due to the 
magnetic nature of the pinning centers §, §. The most 
prominent feature of the vortex lattice pinned by a lat- 
tice of artificial defects is the existence of matching fields 
H n = n$>o/S, where $0 is the flux quantum, and S 
is the area of the primitive cell of the artificial lattice. 
For matching fields, the number of vortices per unit cell 
of the artificial lattice becomes an integer number. As 
expected, such a vortex lattice is pinned much stronger 
by such artificial defects and the film resistivity exhibits 
deep minima (5| at those matching fields. These artifi- 
cial pinning arrays were successfully used to gain insight 
into the macroscopic commensurability effects, but the 
origin and microscopic nature of this phenomena have 



not yet been fully explained. Moreover, the role of the 
self-magnetic field of a ferromagnet in the vicinity of the 
superconductor is not fully understood. 

Therefore, in this letter we consider the interaction be- 
tween a single ferromagnetic disk (FD) and a supercon- 
ducting thin film (SC), within the non-linear GL theory. 
The FD lies on top of the superconductor (xy plane) and 
it is magnetized in the positive z-direction. To avoid the 
proximity effect and exchange of electrons between FD 
and SC we separate them by a thin layer of insulating 
oxide, as is usually the case in experimental conditions. 
Moreover, we consider the magnetic disk to be made of 
a hard magnet whose uniform magnetic moment and in- 
ternal currents are not affected by nearby circulating su- 
percurrents. The ferromagnetic disk is the only source 
of magnetic field and we study in detail how the sys- 
tem is perturbed in the neighborhood of the disk. The 
creation of vortices due to this field, as well as their be- 
havior strongly influence the pinning of additional exter- 
nal flux lines. Furthermore, we find new ordered vor- 
tex/antivortex structures. 

For thin superconductors (d < £, A) it is allowed to 
average the GL equations over the film thickness and 
write them as 

(-zV 2D -^) 2 * = *(f-|4'| 2 ), (I) 

-A 3D l = -^S(z)7 2D , (2) 

where 

?2D = ^ (** Vac* - * V 2D **) - |*| 2 A, (3) 

is the density of superconducting current and A is the 
total vector potential from the FD and supercurrents, 
with boundary condition A = far away from the FD. 
Here the distance is measured in units of the coherence 
length £, the vector potential in ch/2et;, and the mag- 
netic field in H C 2 = cft/2e£ 2 = k\/2H c . The indices 
2D, 3D refer to two- and three-dimensional operators, 
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FIG. 1: The Gibbs free energy as a function of the magnetic moment of the disk with radius Rd/£. = 1 placed on top of the 
superconductor. Top axis shows the flux captured by the superconductor in the positive stray field area of the magnetic disk. 
Insets and figures on the right are the contourplots of the Cooper pair density (red(blue): high(low) density) for a few ground 
state vortex configurations. A giant-vortex is surrounded by antivortices, and the total vorticity equals zero. Only 1/16 of the 
total simulation area is shown. 



respectively. To solve the system of Eqs. we apply 

a finite-difference representation of the order parameter 
and the vector potential on a uniform Cartesian space 
grid (x,y), with typically 512 grid points along the width 
of the rectangular simulation region L x r y \, and use the 
link variable approach [Q, and an iteration procedure 
based on the Gauss-Seidel technique to find 'F. The vec- 
tor potential is then obtained with the fast Fourier trans- 
form technique. The first GL equation is solved with an 
iteration procedure || . The dimensionless Gibbs free en- 
ergy is calculated as F = V^ 1 J{2{A- A )j 2 d - l*) 4 )^, 
where integration is performed over the primitive cell vol- 
ume V, and Aq is the vector potential of the magnetic 
disk. The periodic boundary conditions for A and \& have 
the form j|] 

A{p + h) = A(p) + V Vl (p), (4) 
*(p + 6i) = yexp(2irir h (p)/<P ), (5) 

where bi, i = x,y are the lattice vectors, and rji is the 
gauge potential which cannot be chosen freely but must 
preserve the single valuedness of A and These bound- 
ary conditions mean that A, 4" are invariant under lat- 
tice translations combined with specific gauge transfor- 
mations H . Other quantities, such as the magnetic field, 
the current and the order parameter density are periodic. 

We choose r\ x = Aq x , r\ y = A§ y , with the simulation re- magnetic disk on top 1 1 



gion typically hundred times larger than the magnetic 
disk, i.e. L x >> R c i, which implies the quantization of 
the flux, namely total flux through the superconductor 
equals zero. These boundary conditions result in a pe- 
riodic repetition of not only the superconductor but the 
magnetic disk lattice as well. However, due to the large 
simulation region, these disks are far from each other and 
we are allowed to treat this problem as a single disk on 
top of an infinite superconductor. 



To find the different vortex configurations, which in- 
clude the metastable states, we search for the steady- 
state solutions of Eqs. (|l|) and (Q) starting from different 
randomly generated initial configurations. Then we in- 
crease / decrease slowly the magnetic moment of the mag- 
netic disk m and recalculate each time the exact vortex 
structure. We do this for each vortex configuration in a 
magnetic moment (m) range where the number of vor- 
tices remains constant. By comparing the Gibbs free 
energies of the different vortex configurations we obtain 
the ground state. The results are shown in Fig. [l] for 
the magnetic disk with radius Rd/S, = 1 and thickness 
dd/S, — 0.5, on top of the superconductor with thickness 
<i/£ = 0.5. The magnetic moment is expressed in units 
of mo = cft£/2e = H^S? . Demagnetization effects are 
taken into account in this calculation and the Ginzburg- 
Landau parameter n was chosen to be 1, which approxi- 
mately corresponds to the experimental values found for 
Pb, Nb, or Al films. The values for d, dd and n mentioned 
above will be used throughout this paper. 

From Fig. [l], we notice that with increasing magneti- 
zation of the ferromagnetic disk, the ground state goes 
through different vortex states denoted by successive Ro- 
man numbers. In the case of a finite extend of the su- 
perconducting film, i.e. thin mesoscopic disks, with a 

the Roman numbers corre- 
sponded to the vorticity of the ground state. However, 
in that case, the total flux penetrating the finite size su- 
perconductor was positive, and therefore, the states with 
positive vorticity dominated the free energy diagram. In 
the present case, the magnetic disk is on top of an infinite 
SC film, consequently all flux is captured by the super- 
conductor, and the total penetrating flux equals zero. 
Due to phase conservation, vortices cannot appear indi- 
vidually, but only as vortex-antivortex pairs. The vortex 
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is located under the disk, and the antivortex just out- 
side the positive field region. In Fig. [I], each new vor- 
tex state corresponds to the appearance of a new vortex- 
antivortex pair. Because the magnetic disk is relatively 
small, vortices are strongly confined by the positive mag- 
netic field under the disk, which leads to a giant-vortex 
as the energetically favorable state, while antivortices are 
symmetrically arranged on a ring around the central vor- 
tex (see Fig. Q). The Cooper pair density contourplots 
in Fig. [l] show the states in question for the particular 
value of the magnetic moment of the disk when they be- 
come the ground state, namely at the crossing points 
of the free energy curves. With increasing magnetiza- 
tion of the disk and fixed vortex state, the antivortices 
move further away from the magnetic disk. The config- 
urations shown in Fig. [I] are the result not only of the 
magnetic field of the ferromagnetic disk, but also of the 
interaction between the vortices. While the giant-vortex 
strongly attracts the antivortices, the antivortices repel 
each other. One therefore expects that the vortex con- 
figurations strongly depend on the size of the magnetic 
disk, since it determines the distance between the giant- 
vortex and the antivortices, and the circumference of the 
circle the antivortices are located on. 

Because the vortex configuration is determined by the 
balance between the giant vortex-antivortex attraction 
and the repulsion between antivortices, it is clear that 
the number of antivortices on a ring around the giant- 
vortex cannot increase monotonously. Actually, the for- 
mation of a second, or of even more rings is possible with 
increasing m (see Fig. [j], for the states with N = 9 — 12, 
where N denotes the number of antivortices). A second 
ring is formed when the number of antivortices reaches 
a saturation value n s (or "magic number"), in our case 
this occurs for n s — 6. This, however, does not mean that 
the number of vortices in the first ring remains constant 
with increasing magnetization. For example for N = 12 
it turns out that one more antivortex fits in the first ring. 

Since the antivortices are arranged around the positive 
magnetic field region, the saturation number n s should 
depend on the size of the magnetic disk. Furthermore, for 
larger disks the wider area of the positive field region un- 
der the magnetic disk could allow the giant-vortex under 
the disk to split into individual vortices. Therefore, we 
investigated the influence of the disk size on the vortex 
structure of the superconducting film. We enlarged the 
magnetic disk from Rd/S, — 0.5 to 5 keeping the thick- 
ness dd and the parameters of the SC fixed. In doing so, 
we obtain the equilibrium vortex phase diagram, shown 
in Fig. ||. Solid lines correspond to transitions between 
the successive N states, and the thick dashed lines give 
the number of antivortices needed for the formation of 
an additional ring of antivortices. In the white (shaded) 
region the giant- (multi-) vortex state is found under the 
magnetic disk. 

Notice that with enlarging the magnetic disk the num- 
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FIG. 2: Phase diagram: the relation between the radius of 
the disk, its magnetic moment, the number of antivortices 
(N), the number of antivortex rings (Roman numbers) and 
the vortex state under the magnetic disk (shaded area: mul- 
tivortex state). 

ber of antivortices needed before a new ring is formed 
increases, as expected. But, on the other hand a larger 
magnetic moment is needed for the appearance of a new 
vortex-antivortex pair in the ground state. However, the 
volume of the disk also increases, which diminishes the 
effects of the magnetic moment increase. Therefore, we 
calculated the flux <f> + through the SC, due to the pos- 
itive part of the stray field in the z = plane, namely 
in the center of the superconductor and up to the ra- 
dius i?o— Rd where the magnetic field changes its polarity. 
The top axis in Fig. [l] shows the value of this flux. Such 
a calculation tells us that if we keep the magnetic mo- 
ment constant and increase the radius (volume) of the 
disk, the flux <I> + decreases, making the positive (and 
negative) flux through the SC approximately constant 
along the transition line between the successive states in 
Fig. |2|. Obviously, a larger flux is needed to create the 
first vortex-antivortex pair, i.e. A$ + /$o = 2.112, due 
to the high stability of the Meissner state and the asym- 
metry of the N = 1 state (see inset of Fig. [j]). Further 
increase of the flux $ + decreases A$ + and for larger 
N it approaches the value A$+/$ = 1.073. Notice 
that the flux <I> + is not exactly quantized in units of $o 
which is a mesoscopic effect. The quantization condition 
$ = § c A ■ dl = L$o cannot be used because it is not pos- 
sible to construct a contour C around the positive stray 
field region where the current is zero. 

With increasing radius of the disk, the giant-vortex 
splits into individual vortices, as shown in Fig. ||(a). 
Figure ||(b) shows the phase of the superconducting or- 
der parameter, from which the location of the vortices 
and antivortices can be easily deduced. With increas- 
ing magnetic field, these vortices move closer together, 
resulting ultimately in the formation of a giant-vortex. 
Figs. |^(c,d) show the current distribution in the SC: the 
supercurrents around each vortex and antivortex flow in 
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FIG. 3: Splitting of the giant- vortex under the disk into a mul- 
tivortex: a) Cooper pair density, b) superconducting phase 
contourplot (zero - blue, 2ir - red), for N = 3, Rd/£ = 4.0, 
and m/mo = 18.5. (c) is a vectorplot of the corresponding 
supercurrents and (d) is an enlargement of the central area. 

opposite directions, which are determined by the direc- 
tion of the change of the phase of the order parameter. 
However, the current around the central vortices (i.e. at 
the edge of the magnetic disk) has the same direction 
as the one around the antivortices, although the phase 
change implies otherwise. The reason is that the direc- 
tion of this supercurrent is governed by the vector poten- 
tial of the magnetic field of the ferromagnetic disk which 
is maximal near the edge of the disk. The value of the 
total current will be maximal in the same region. This 
is also the reason why with increasing magnetic moment 
of the FD new vortices nucleate exactly under the edge 
of the disk, i.e. when the maximal current reaches the 
value of the GL current. 

To analyze the generation of new vortices, we start 
from one specific initial condition, for example, the 
ground state configuration for N = 3. Then we increase 
slowly the magnetic moment of the magnetic disk m and 
recalculate the vortex structure within 10 5 iteration steps 
for each value of m. This way we can find the magnetic 
moment for which the initial state changes from N = 3 
to the N — 4 configuration and we observe the nucle- 
ation of a vortex-antivortex pair by plotting the inter- 
mediate results during the different iteration steps. The 
result for the III— ► IV transition is shown in Fig. ||. The 
Cooper pair density plot Fig. |J(a) shows only a distorted 
central vortex, but the superconducting phase plot (b) 
reveals the real vortex structure. Namely, in addition 
to the configuration similar to the state III, one can no- 
tice a vortex-antivortex pair nucleating under the edge of 
the dot. Further iterations in our calculation show that 
the vortex from this pair will move towards the central 
giant-vortex and the antivortex moves to the periphery, 




FIG. 4: Generation of the vortex-antivortex pair with increas- 
ing magnetic moment of the disk: a) Cooper pair density, 
and b) superconducting phase contourplot, for N = 3(4) and 
Rd/S, = 2.0. The thick ring shows the edge of the FD. The 
upper inset is an enlargement of the phase of the central area 
which more clearly shows the vortex-antivortex pair. 



causing the rearrangement of the shell of antivortices into 
the configuration IV. 

The predicted new vortex configurations can be ob- 
served experimentally by using e.g. scanning probe tech- 
niques like Hall and Magnetic Force Microscopy. Further- 
more, these vortex structures will influence the pinning 
properties and the superconducting phase diagram |1T|] . 
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